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The single-particle spectral functions in asymmetric nuclear matter are computed using the ladder 
approximation within the theory of finite temperature Green's functions. The internal energy and 
the momentum distributions of protons and neutrons are studied as a function of the density and 
the asymmetry of the system. The proton states are more strongly depleted when the asymmetry 
increases while the occupation of the neutron states is enhanced as compared to the symmetric 
case. The self-consistent Green's function approach leads to slightly smaller energies as compared 
to the Brueckner Hartree Fock approach. This effect increases with density and thereby modifies 
the saturation density and leads to smaller symmetry energies. 

PACS numbers: 21.65.+f, 21.30.Fe 



I. INTRODUCTION 

The equation of state (EOS) of asymmetric nuclear 
matter is a necessary ingredient in the description of as- 
trophysical environments of interest such as supernova 
explosions or the structure of neutron stars Q . Actually, 
the study of asymmetric nuclear matter is also relevant 
to understand stable nuclei because they themselves are 
asymmetric nuclear systems with a different number of 
protons and neutrons. The recent availability of data 
concerning nuclei far form the stability valley has also 
renovated the interest in the study of asymmetric nu- 
clear matter as a first step in the microscopic study of 
these nuclei. 

The evaluation of this EOS starting from realistic mod- 
els of the nucleon-nucleon (NN) interaction is still one 
of the challenging open problems in nuclear physics. In 
fact, the presence of strong-short range and tensor com- 
ponents in the realistic NN interactions, which are re- 
quired to fit the NN scattering data, are the origin of 
the corresponding correlations in the nuclear wave func- 
tion. The study of these correlations and their influence 
on different observables has recently made an important 
progress not only from the theoretical side but also from 
the experimental point of view 0- C3- DD] - In the later case, 
the analysis of (e, e'p) reactions on 208 Pb covering a wide 
range of missing energies lead to the conclusion that the 
occupation numbers for the deeply bound proton states 
are depleted by about 15-20 % This depletion can be 
identified with the corresponding depletion of hole states 
in nuclear matter with momenta well below the Fermi 
momentum 00, 13- 

Several theoretical tools have been developed and ap- 
plied to describe these correlations in nuclear systems. 
These include the Brueckner hole-line expansion [9| and 
also variational approaches using correlated basis func- 
tions nsm 03. 

Recently, an enormous progress has been achieved 



in using the self-consistent evaluation of Green's func- 
tion Il3l U4[ to solve the nuclear many-body problem 
H H ll O [H H E 13 . This method gives direct 
access to the single particle spectral function, i.e. to the 
single particle properties and, in particular, to the occu- 
pation numbers. At the same time, the modification of 
the single particle properties affects the effective interac- 
tions between nucleons in the medium and both things 
should be determined in a self-consistent way. 

Most of the microscopic calculations have been ad- 
dressed to study symmetric nuclear matter (SNM) and 
pure neutron matter (PNM) [!llll7| . The study of asym- 
metric nuclear matter is technically more involved and 
only few Brueckner-Hartree Fock calculations are avail- 
able [2TI I22I |23| . In most of the cases, one assumes a 
quadratic dependence of the energy per particle 
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in terms of the asymmetry parameter, a = (N — Z)/A 
and the symmetry energy a s (p). In this way, the calcu- 
lation of the energy of SNM and PNM allows to deter- 
mine the symmetry energy a s (p) and using the previous 
equation one can estimate the energy for any asymmetry. 
That this quadratic expression is a good approximation 
has been directly confirmed in Brueckner-Hartree-Fock 
calculations of asymmetric nuclear matter [2lH22| . In the 
case of the variational approach, calculations for asym- 
metric matter with the same accuracy that the ones per- 
formed for SNM or PNM are not yet available. However, 
besides the energy per particle which is governed mainly 
by the symmetry energy, there are other observables, as 
the momentum distributions, which do not necessary fol- 
low a quadratic dependence. 

In this paper we want to perform a calculation of asym- 
metric nuclear matter in the framework of SCGF theory. 
The calculation is performed at finite temperature. A 
temperature of 5 MeV has been chosen because it is small 
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enough to allow for conclusions for the T = case. On 
the other hand this temperature is large enough to al- 
low for a smooth numerical representation of the spectral 
functions and to avoid the pos sibility of proton-neutron 
pairing instabilities |2 lL l25l l26| . In the next section, we 
briefly describe some specific features of the SCGF for- 
malism for asymmetric nuclear matter. The discussion 
of the results and a systematic comparison with results 
of the Brueckner-Hartree-Fock approach is presented in 
Section III. The main conclusions are summarized in the 
last section. 



II. FORMALISM 

A key quantity in the theory of Green's functions - 
that allows to compute all single-particle observables as 
well as the internal energy of the system — is the spectral 
function A{k,ui). It can be obtained as a solution of 
Dysons's equation, 
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where £ T (fc, w+irj) is the retarded self energy of a nucleon 
with isospin projection r, which can either be a proton 
(t = +| or p) or a neutron (r = — \ or n). 

In the ladder or T matrix approximation, the self en- 
ergy contains an energy independent Hartree-Fock part 
and a complex dispersive contribution that accounts for 
correlations between the particles. The Hartree-Fock 
contribution to the ladder self energy of a nucleon r in- 
volves an explicit sum over the isospin projection t' of 
the internal particle, 
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with V^fT the nuclear two-body potential in a partial 
wave representation. The relative momentum between 
the interacting particles is given by q = 4(k — k') and 
the Clebsch-Gordan coefficient for the isospin quantum 
numbers is denoted by Cf T f T , . The momentum distri- 
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must be derived from the non-trivial spectral function. 
The Fermi-Dirac function for nucleons with isospin pro- 
jection r is denoted by f T (u>) = {exp[/3(w — /i T )] + 
where (3 stands for the inverse temperature T _1 . For a 
given total density p, the partial fraction x T = p T / P of 



the respective particle species is given by 
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where 7 = 2 is the spin degeneracy factor of the system. 
The partial fractions of protons and neutrons add up to 
one, x p + x n = 1 and the asymmetry is given by a — 
x n — x p . Considering a fixed composition, Eq. (jjjj can 
be used to fix the chemical potential fi T . 

The generalized expression for the imaginary part of 
the self energy in the T matrix approximation for the 
case of asymmetric nuclear matter reads: 
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Here, b T y(Q) = {exp[/3(fi 
Bose distribution function 
is given by P — i | k + k' | . 
is reported, e.g., in Ref. 
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fi T — fi T ')] — is the 
The total pair momentum 
A dispersion relation that 
determines the real part 



of the dispersive contribution to the ladder self energy. 
The T matrix elements contain the re-summation of the 
ladder diagrams to all orders and can be obtained as the 
solution of a scattering-type integral equation, 
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Since in ANM, the in-medium propagation of a neutron 
is different from that of a proton, the T matrix elements 
depend upon the third component of the isospin of the 
propagating pair, tjit = t + r' , via the non- interacting 
two-particle propagator g\} %T , 



g^ i+T2 (k 1 ,k 2 ,n + i V )= (8) 
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To circumvent the coupling between partial waves of dif- 
ferent total angular momenta, g^ T enters Eq. Q in an 
angle-averaged form that is indicated by the bar. 

Eqs. J3J |[SJ must be solved self-consistently for a given 
temperature /3 _1 , a given density p and a given proton 
fraction x p . Partial waves up to J = 9 were included for 
the calculation of the energy-independent part of the self 
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energy. The dispersive part includes only partial waves 
up to J = 2. A more comprehensive description of the 
numerical details is given in Ref 0|- The numerical 
routine that was developed for ANM was tested in the 
following way: the results for SNM and PNM are re- 
covered if one chooses x p — 0.5 and x p — 0, respec- 
tively. Furthermore, we have checked that charge sym- 
metry is fulfilled, which can be expressed by the condition 
T, p (k,uj,Xp) = E n (fc,w, 1 - x p ). 

Once the self-consistent solution for the spectral function 
has been obtained, the internal energy per particle in the 
system is given by 

T — oo 

(9) 

The BHF approach is obtained from the previous for- 
mulation by assuming that the single particle spectral 
functions are characterized by only one energy having 
the full strength accumulated in this energy, A T {k, w) = 
S(uj - e^ HF (k)), being e^ HF (k) the BHF single particle 
energy. In addition, in the two-body propagator (Eq.lJSJ) 
one only considers the propagation of particle states. 
Notice also that the BHF self-energy does not include 
the contribution of the Bose distribution that appears in 
Eq.©. 



III. RESULTS AND DISCUSSION 

All the results discussed in this paper have been com- 
puted for the charge-dependent Bonn (CDBONN) po- 
tential, defined in [27j which is non-local and exhibits a 
softer tensor component compared to other realistic po- 
tentials as the Argonne V18 or Reid93 [2{j, which are 
local. Since we want to concentrate on the asymmetry 
dependence, we will only consider one temperature T = 5 
MeV, low enough for the conclusions on the asymmetry 
dependence to be valid at T = MeV and high enough 
to avoid the instabilities associated to the neutron-proton 
pairing. 

The binding energies of asymmetric nuclear matter cal- 
culated in the BHF approximation and in the framework 
of the SCGF theory are reported in Fig. as a func- 
tion of the square of the asymmetry parameter a 2 . The 
plots correspond to two densities, the empirical satura- 
tion density of symmetric nuclear matter, p = 0.16 fm , 
and twice this density. Notice, however, that BHF calcu- 
lations using the CDBONN potential yield a saturation 
point at higher densities. 

The first thing to realize is the linear dependence of 
the energy in terms of a 2 , for both types of calculations 
(BHF and SCGF) in the full range of variation of the 
asymmetry from SNM to PNM. In the case of BHF, this 
fact has already been considered in the literature [2l| 
and provides the justification to use PNM and SNM to 
determine the coefficients in Eq. which turns out to 




0.2 0.4 0.6 0.8 



Figure 1: Dependence of E/A on a 2 computed in the SCGF 
approach (full lines) and in the BHF approach (dashed lines) 
at two densities, p = 0.16 fm -3 (black lines) and p = 0.32 
fm~ 3 (red lines). 



be a very good approximation. 

The propagation of holes in the SCGF and the con- 
sideration of the spectral functions in the intermediate 
states of the ladder equation result in a repulsive ef- 
fect with respect to the continuous choice BHF calcu- 
lation. This repulsive effect increases with density, being 
1.8 MeV and 3.2 MeV in SNM matter at p = 0.16 fm" 3 
and p — 0.32 fm -3 , respectively. 

However, for a given density, the difference between 
the SCGF and the BHF calculation does not depend so 
much on the asymmetry, being slightly larger in SNM 
than in PNM. This leads to a small decrease of the sym- 
metry energy of the SCGF calculation as compared to 
the BHF one. At p = 0.16 fm" 3 , a s = 30.0 MeV in 
the BHF approximation while the SCGF schemes provide 
a s = 28.6 MeV. This is in qualitative agreement with a 
recent comparison of both approaches \u[ , in which the 
symmetry energy has been deduced from SNM and PNM 
calculations. The reduction of the symmetry obtained 
in the present analysis, however, is weaker than the de- 
crease by around 4 MeV reported in Il7j . In this refer- 
ence, however, the Reid93 potential, which has a stronger 
tensor component and also stronger short range correla- 
tions than the CDBONN, was employed. Furthermore it 
should be noted that the spectral functions employed in 
|17| are described in terms of three 5 functions, whereas 
a continous description has been used in the present ap- 
proach. In addition, the calculation of was performed 
at T = whereas we consider T = 5 MeV. We do not 
believe, however, that the temperature can be the origin 
of this discrepancy. In fact, the effect of the finite tem- 
perature is a slight increase of the symmetry energy, but 
this will also be true in a similar amount for the BHF ap- 
proach. Instead this discrepancy demonstrates that the 
loss of energy in SCGF as compared to BHF is a delicate 
balance between an increase of kinetic energy and more 
attractive potential energy. In fact, the perturbative in- 
clusion of the hole-hole scattering terms of |2^| even leads 
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Figure 2: Dependence of neutron (upper curves) and proton 
(lower curves) chemical potentials p, on the asymmetry a cal- 
culated in the SCGF approach (solid lines) and on the BHF 
approach (dashed curves) at two given densities, p = 0.16 
fm -3 (black lines) and p — 0.32 fhi _3 (red lines). 



to an increase of the symmetry energy. 

The second point that we want to consider is the de- 
pendence of the chemical potential of protons and neu- 
trons on the asymmetry. Fig. (0) shows the chemical 
potential of neutrons (upper curves) and protons (lower 
curves) calculated in the SCGF approach (solid lines) 
and in the BHF approach (dashed curves). The chem- 
ical potential of protons and neutrons coincide at the 
symmetric case (a = 0), as we have not considered the 
charge symmetry breaking terms contained in the CD- 
BONN potential. When the asymmetry increases, the 
neutron chemical potential increases and becomes posi- 
tive, while the one of the protons becomes more and more 
attractive. The dependence on a 2 is not linear anymore. 
One should keep in mind that in the case of SCGF, the 
chemical potential that one obtains from the normaliza- 
tion condition of the partial density should coincide with 
the one obtained from the free energy using the thermo- 
dynamic relation p, T — F/A + p ®^ p A ■ On the contrary, 
it is well known that in the BHF approach the chemical 
potential derived from the normalization condition sub- 
stantially differs from the one obtained by the thermo- 
dynamic relation. Actually, this is already true at T = 
MeV where one needs to incorporate the rearrangement 
terms in the self-energy in order to recover the relation 
fi = e(kf) (where e(fcp) is the single-particle energy at 
the Fermi surface). 

It is important to note that the differences between the 
BHF and the SCGF chemical potentials are much larger 
(of about ~ 15 MeV in the whole range of asymmetries) 
than the differences in the energy per particle discussed in 
Fig. , thus indicating that the role of correlations has a 
larger influence in this observable than in the energy per 
particle and also giving an idea of the magnitude of the 
rearrangement term in such calculations. In particular 
we want to emphasize that the difference between the 
chemical potentials for neutrons and protons is larger in 
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Figure 3: k — MeV proton spectral function for different 
proton fractions. The calculations were performed at p = 0.16 
fm _3 and T = 5 MeV. 



the case of BHF as compared to SCGF. This implies that 
the SCGF tends to predict /3-equilibrium with a smaller 
proton fraction as derived from BHF calculations. 

The next point we shall address is the discussion of 
the single-particle spectral functions. Fig. @ shows the 
proton spectral functions at p = 0.16 fm - for the mo- 
mentum k = MeV at different proton fractions as a 
function of the energy measured with respect to the pro- 
ton chemical potential corresponding to each fraction. 
As the asymmetry increases, the amount of protons de- 
creases and the Fermi-momentum for the protons gets 
closer to k — MeV. As a consequence, the coupling 
to two-hole one-particle configurations with this momen- 
tum of k = is reduced, the quasi-particle peak obtains 
a smaller width and the peak gets higher. The spectral 
function at positive energies, however, is larger with in- 
creasing asymmetry. This can be explained by the fact 
that correlation effects mainly originate from the ten- 
sor force, which is more important in the proton-neutron 
than in the neutron-neutron or proton-proton interac- 
tion. Note that the most important partial wave for these 
tensor correlations is the 3 Si — 3 Di channel, which is rel- 
evant for nucleon pairs with isospin zero only. Since the 
density of neutrons increases with the asymmetry pa- 
rameter, the protons display stronger correlation effects 
at these larger asymmetries leading to an enhancement 
for the spectral function for k < kp at energies to > p. 

The corresponding plot for the neutron spectral func- 
tions is presented in Fig. JIJ|. Again we consider the 
momentum k = 0, the same density p — 0.16 fm -3 and 
asymmetries. Looking at the spectral functions at en- 
ergies above the chemical potential for the neutrons, we 
can see that the correlations which are responsible for 
the spectral function in this regime are reduced with the 
asymmetry, i.e. with the density of protons. This can 
again be understood from the dominance of the proton- 
neutron interaction leading to these correlations. The 
width of the neutron spectral functions seems not to be 
affected very much by an increase of the asymmetry for 




Figure 4: k = MeV neutron spectral function for different 
proton fractions. The calculations were performed at p = 0.16 
fm -3 and T = 5 MeV. 
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Figure 5: Dependence of the k = MeV on-shell width T(0) 
on the proton fraction for neutrons (full lines) and protons 
(dashed lines) at two given densities, p = 0.16 fm -3 (black 
lines) and p — 0.32 fm _3 (red lines). 



the considered values of a. This is a result of the fact 
that the damping of the strong isospin zero correlations 
at larger values of a is counterbalanced by an increasing 
phase space for the neutron- neutron configurations. Sim- 
ilar observations have also been made in the perturbative 
calculations of asymmetric matter in |23| . 

In order to explore a bit more the width and height of 
the quasi-particle peak in the spectral function and its 
dependence on the asymmetry we show in Fig. (jSJ) the 
width of the spectral function given by 2 | 7ro£(fc, e qp ) \, 
where e qp is the quasi-particle energy. Obviously, for the 
symmetric case both the widths for protons and neu- 
trons coincide. When the asymmetry increases, the width 
of the protons decreases monotonically, indicating that 
k = is getting closer and closer to the proton Fermi 
surface. The situation for the neutrons is a little bit 
different. For low asymmetries, the width initially in- 
creases with the asymmetry, reaches a maximum and 
then it decreases again - being a little lower for pure 
neutron matter than for the symmetric case. This in- 



Figure 6: Dependence of the n(0) on the proton fraction for 
neutrons (full lines) and protons (dashed lines) at two given 
densities, p = 0.16 fm _3 (black lines) and p = 0.32 fm _3 (red 
lines). 



dicates that there is an asymmetry for which the corre- 
lation effects on the width of the neutron quasi particle 
peak reaches a maximum. As mentioned above, there is 
a competition between the correlations originated by the 
propagation of neutron-neutron holes, which should in- 
crease when the asymmetry gets larger as the associated 
phase space also grows, and the correlations originating 
from the interaction of the neutrons with the protons. 
It is clear that these last correlations get less important 
for the neutrons with increasing asymmetry. Finally, for 
neutron matter, the decrease of the number of protons 
and therefore the suppression of the tensor components of 
the neutron-proton interaction dominates and the width 
for neutron matter is slightly smaller than for symmetric 
nuclear matter. 

Another observable which has been object of discus- 
sion during the last years, due to its relevance in the 
analysis of the (e, e'p) reaction is the momentum distri- 
bution, which also provides a clear measurement of the 
effects of correlations. In Fig. iJBJ, we show the occu- 
pation of the zero momentum state as a function of the 
asymmetry. In the symmetric case, we can observe an 
unexpected behavior on the occupation as a function of 
density: at the larger density there is a higher occupation 
while at the lower one the occupation is smaller, a behav- 
ior that has been pointed out previously in the literature 
For a given density, the proton depletion increases 
with the asymmetry, indicating the importance of the 
neutron-proton correlations. On the contrary, the deple- 
tion of the neutrons gets smaller, i.e neutron matter is a 
less correlated system. Also it is worth to notice that for 
pure neutron matter one recovers the expected behavior 
of having a higher occupation for smaller density. 

Finally, in Fig. J7J, we show the momentum distribu- 
tion at p = 0.16 fm~ 3 and an asymmetry a = 0.2, which 
is characteristic of heavy nuclei, such as 208 Pb. Included 
in the figure are also the BHF momentum distributions. 
The first thing to notice is that the BHF momentum dis- 
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Figure 7: Momentum distributions for neutrons and protons 
for the SCGF approach (full lines) and the BHF approach 
(dashed lines) at a density p — 0.16 fm -3 and asymmetry 
a = 0.2. 

tributions do not contain correlation effects and are very 
similar to a normal thermal Fermi distribution. The mo- 
mentum distributions obtained in the framework of the 
SCGF contain, besides the thermal effects, important 
short-range and tensor correlations which are reflected 
in the depletion of the occupation at low momentum and 
in a larger occupation than the BHF momentum distribu- 
tion at large momenta. Notice also that the proton mo- 
mentum distribution is more depleted than the neutron 
momentum distribution. This is in agreement with all the 
previous discussions, reflecting the fact that the protons 
(i.e. the less abundant particle in the realistic asymmet- 
ric conditions) are more affected by correlations, mainly 
due to its interactions with neutrons. This may be im- 
portant in the interpretation of the momentum distribu- 
tions obtained in the (e, e'p) experiments. Some of the 
analysis have been conducted in the framework of a lo- 
cal density approximation starting from results obtained 
in SNM, where the momentum distribution of protons 
and neutrons are identical. This different behavior of 
the momentum distributions of protons and neutrons is 
in contrast with the very recent calculations by Bozek 
[20| obtained also in the SCGF framework and with the 
CDBONN potential in which no noticeable difference be- 
tween the momentum distributions of neutrons and pro- 
tons was found. 



IV. SUMMARY AND CONCLUSIONS 

The techniques to evaluate the single-particle Green's 
function in a self-consistent T- matrix approach (SCGF), 



that has recently been developed for nuclear matter [16( , 
have been extended and applied to asymmetric matter. 
Actual calculations have been performed using the re- 
alistic CDBONN interaction for asymmetric matter at 
two densities (the saturation density of symmetric nu- 
clear matter and twice this density) , various asymmetries 
and for a temperature of T = 5 MeV. This temperature 
is low enough to allow for conclusions on the T — limit 
and high enough to avoid the instabilities associated to 
the neutron-proton pairing. 

The inclusion of the hole-hole ladders and the self- 
consistent treatment of the Green's Function in the 
SCGF approach leads to a small reduction of the binding- 
energy per nucleon as compared to the BHF approxima- 
tion. This effect increases with density and is slightly 
weaker for pure neutron matter as compared to symmet- 
ric nuclear matter which leads to a small reduction in the 
symmetry energy. 

Larger effects are observed for the single-particle prop- 
erties like the chemical potential. In particular we ob- 
serve in neutron rich matter a reduction in the difference 
between the chemical potential for protons and neutrons, 
which would correspond to smaller proton fractions in /3 
stable matter than those predicted by a BHF-type equa- 
tion of state. 

The SCGF calculation also yields detailed informations 
on the single-particle spectral functions and momentum 
distributions for protons and neutrons. We observe a 
depletion of the proton occupancies for momenta below 
the Fermi momentum, which increases significantly with 
the neutron fraction. This can be explained by the strong 
correlations induced from proton neutron interaction. 
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